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The chiral phase transition in the conventional random matrix model is the second order in the 
chiral limit, irrespective of the number of flavors Nf, because it lacks the UA(l)-breaking deter- 
minant interaction term. Furthermore, it predicts an unphysical value of zero for the topological 
susceptibility at finite temperatures. We propose a new chiral random matrix model which re- 
solves these difficulties by incorporating the determinant interaction term within the instanton gas 
picture. The model produces a second-order transition for Nf ^2 and a first-order transition for 
Nf ^3, and recovers a physical temperature dependence of the topological susceptibility. 
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1. Introduction 

One of the prominent features of non-perturbative QCD is spontaneous breaking of chiral 
symmetry in the light quark sector. According to the Banks-Casher relation [[I]], the order parameter 
{qq) is related to the Dirac spectral density at zero eigenvalue in the thermodynamic limit. Within 
the so-called e regime characterized by the system size such that niji <^ l/L <^ 1 GeV, the 
QCD partition function is dominated by the constant pion field configurations. In this regime, 
universal properties of the Dirac eigenvalue distribution near zero are legitimately analyzed in a 
chiral random matrix (ChRM) theory [Q], where the kinetic term of the Dirac operator is discarded 
and the complexity of the gauge field dynamics is represented by treating the Dirac operator as 
a random matrix of constant modes. The matrix size of the constant modes 2N is considered to 
be proportional to the system volume. By taking a large volume limit N ^ o°, away from the £ 
regime, one can study thermodynamics of the ChRM theory as a schematic mean-field model for 
QCD. One finds the ground state of the model in the chirally broken phase and can explore the 
model phase diagram by introducing the temperature t [^] and the quark chemical potential jJ. [Q]. 

There are two drawbacks, however, in the ChRM model, concerning the Ua(1) symmetry. 
Firstly, the UA(l)-breaking determinant term is missing in the effective action of the model, which 
consequently predicts a second-order phase transition at finite temperature for any number of the 
quark flavors Nf. Secondly, the topological susceptibility is suppressed unphysically to vanish at 
finite temperatures. 

In an earlier work [Q], the appropriate form of the ChRM model with a UA(l)-breaking term 
is speculated from a quark model with the determinant interaction in O-i-l dimensions. In addition 
to the constant modes in the conventional ChRM model, new constant modes are introduced so as 
to reproduce the determinant interaction. These new modes are considered to be associated with 
instantons and called topological zero modes. However the effective potential of the starting quark 
model is unbounded from below for Nf = 3, and therefore no physical ground state exists. 

In this paper we propose a new ChRM model [^] changing the distribution of the topological 
zero modes, which results in the determinant term appearing under a logarithm. Then the effective 
potential becomes stable for any Nf, and describes a second-order phase transition for Nf = 2 while 
a first-order transition for Nf = 3. Moreover, the topological susceptibility shows physical behavior 
as a function of temperature in our model. 



2. Chiral random matrix model at finite temperature 

In analogy with the QCD partition function, the ChRM model with Nf flavors of mass m f is 
defined as an average of the quark determinants |^ 

Zv = / dWe-^^'''^'^ Yl det(D + mj) , (2.1) 

where the Dirac operator has been replaced with an anti-Hermite matrix of constant modes 

D=( ' + 

I /WU/flw-|v|/2 i 
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with an (A'^ + v/2) x {N — v/2) complex matrix W. The Gaussian random distribution of W is 
a simple realization of the complexity of the gluon dynamics. Chiral symmetry is retained as a 
fact {D, 75} = with 75 = diag(l;v+v/2i ~^n-v/2)- The temperature effect has been introduced in 
the Dirac operator with a constant t, which may be interpreted as the lowest Matsubara frequency 
nT. One can easily show that this matrix D has |v| exact zero eigenvalues, which are interpreted 
as the zero modes accompanied by the topological charge V. The complete partition function is 
obtained after the sum over v weighted by the quenched topological susceptibility T of the pure 
gluon dynamics: 

Z0= ^ 2P'v)-e"'^Zv , (2.3) 

v=-2N 

where the 6 angle has been introduced. 

After rewriting the determinant with the fermion variables and doing the Gaussian integral 
of W in Zy, we find a four-fermion vertex interaction, which can be unfolded by introducing an 
Nf X Nf auxiliary variable S ~ q^qR- We perform the fermion integral to obtain an expression for 
Zy in terms of S, 



Z. = /^5.--M..^)det((5 + ^)(5t + ^)+.2)^-^x|'^^^' + -^)^ ^'-'^ , (2.4) 
J ;v ; ; \det(5l' + ^1")-^ (v < 0) 

where ^ = diag(mi , . . . ,mNf). 

In the thermodynamic limit we set v = and evaluate Eq. (2.4) with the saddle point equation. 
For Soc 1^^ and ^ = 0, the Nf dependence is factored out in the saddle point equation, implying 
a second order phase transition for any number of Nf. Furthermore, Zy is non-analytic in v as 
the integrand contains a term with | v | , which causes the unphysical suppression of the topological 
susceptibihty at finite temperatures [^. Note that |v| disappears when t = 0. 



3. Chiral random matrix model with determinant interaction 

Near-zero modes and topological zero modes - Let us first recall the instanton gas picture. An 
isolated instanton is a localized object accompanying a right-handed exact fermion zero mode. In 
a dilute system of N+ instantons and N- anti-ins tantons, we expect A'^+ right-handed and left- 
handed zero modes even at finite temperatures. In an effective theory at long distances, effects of 
the instantons should be integrated out, which will result in UA(l)-breaking effective interactions. 
The fundamental assumption in our model is the classification of the constant modes into the near- 
zero modes and the topological zero modes ^ ^. We deal with the 2N near-zero modes appearing 
in the conventional models and additionally the N+ + A'^- topological zero modes which we regard 
as the modes accompanied by the instantons. In our model, N± fluctuate according to the instanton 
distribution, and the number of the exact zero modes is given by v = A'^^ — A'^-. Eventually we sum 
over N+ and N- assuming a certain distribution with the mean value of &{N). 

We write down a Gaussian ChRM model for definite numbers of zero modes as [|5|, |6|] 

< ^ = / JBJX JFe-^^'"-(^'+^^' +^^'+^^') n det(£» + nif) (3. 1) 
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with 



D 



( /A + jYIat iX\ 
iA^ ^it\t^ iY 
/y^ iB 

\ iX^ iB'f J 



(3.2) 



where the matrix A corresponds to the near-zero modes of the conventional model and the A'^^ x N- 
matrix B represents the topological zero modes. The matrices X and Y induce mixing among these 
modes. Note that the temperature term t is introduced only for the near-zero modes. Following the 
same steps as in Sec. 0, we find the sigma model representation for this ChRM model, which is 



analytic in N± as well as A'^ unlike Eq. (2.4): 



7^ 



dS e 



-NI.-trS''S 



det [{S + ^) (51" + + ^ det(5 + det(5^ + 



(3.3) 



Distribution of the topological zero modes - The complete partition function is obtained after 
summing Z^^ over the instanton numbers A'^^ and A'^- . Here we simply assume independent 
distributions P{N±) for N+ and A'^-, i.e., 



Ze= e'''PiN+)PiN^)Zl,^_ = fdSe-'^''^'''-"^''K 



(3.4) 



Let us first consider P{N±) in a dilute instanton gas picture. For a one-instanton configuration, one 
may assign a weight K compared with a no-instanton configuration, and multiply a factor N °<^V 
taking into account the integration over the instanton location. For a configuration with A^+(_) 
(anti-)instantons, we then have a Poisson distribution 



Ppo(A^±) 



1 

A^ 



[kN) 



N± 



(3.5) 



where the factorial A'^_|_(_)! appears as the symmetry factor. The summation with Ppo[N±) in 



Eq. (3.4) results in the exponentiation of the determinant term 



Po 



■ ' r^tr^^t _ 1 In det [{S + ^) (5^ + + - i Kie''^ det(5 + 



(3.6) 



This determinant term is commonly incorporated in effective models as the Ua(1) anomaly term. 
However this potential is unbound for Nf = 3 in the ChRM model because the term det(5 + ^) ~ 
<p^ for large S = (pl^f dominates over the other terms in Q.po. 

It should be noticed here that the fermion coupling distorts the N± distribution itself. With 



including the determinant term of the topological zero modes in Eq. ( |3.3[ ), the effective distribution 
for A'^^ reads 



Ppo{N+) 

with d = I det (5 + ^) \ , and similarly for A^_ 
indefinitely with increasing J ~ 0^^^ as {N±) 



1 



(KNd) 



(3.7) 



This means that the average value of N± increases 
= KNd. However, the possibility of infinitely many 
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Figure 1: Chiral condensate ^ as a function of t and m for Nf = 2 (left) and 3 (right). 



degrees of freedom N± within a finite volume is inadequate for a regularized low-energy effective 
theory. We need a cutoff for N±. 

Here we set explicitly a maximum value of i&{N) forN±. We split the finite space-time volume 
into jN cells with / being a constant of i^(l), and assign a probability p for a cell to be occupied by 
a single (anti-)instanton and (1 — p) for a cell unoccupied. This assumption results in the binomial 
distributions for A^±: 

p^^i^-py"^"^- (3-8) 

For a small p and a large jN, the binomial distribution P{N±) is accurately approximated with the 
Poisson distribution with the mean jNp. But it cannot be a good approximation for a large p. The 
binomial distribution provides a stringent upper bound jN for the number of modes N±, in contrast 
to the Poisson distribution. The corresponding effective potential for S is found to be 



a{S;t,m, e) =^I?trSS'' - ^ Indet [{S + ^) {S^ + ^^) + 1^] 



In (e'^ad&t{S + J^) + 1 ) + ln (e-'^aA&t{S'^ + ^^) + 1 



(3.9) 



with a = p/{\— p). 

The variance of the topological charge v = A'^^ — A'^- for the binomial distribution is computed 
as INt = 2Nyp{ \ — p), where t is the quenched topological susceptibility. In the presence of the 
fermion coupling, this susceptibility will be replaced with 

4. Ground state and fluctuations 

In this section we shall study ground state properties of the system with equal mass, ^ = 
ttiInj, for simplicity. Setting S = (plNf with real (p and with = 0, we obtain a simple form of the 
grand potential: 

D.{(l);t,m) =^NfI?^^ - ^Nfln[{^ +mf + t^] - ^Y^n\a{^ +mff + \\^ . (4.1) 
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Nf=2 Ni = i 
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Figure 2: Temperature dependence of the mesonic masses in the flavor-singlet scalar (sO) and pseudo-scalar 
(psO) channels and in the flavor-nonsinglet scalar (s) and pseudo-scalar (ps) channels for m = (thick lines) 
and m ^ (thin lines). We set m = 0.1 (m^. = 0.0265) as nonzero quark mass for A'/ = 2 (3). 



The factor Nf cannot be factored out in the potential Q. because of the anomaly term here. 

In the thermodynamic limit we calculate the quark condensate {qq) oc for £ = 1, a = 0.3 and 
7 = 2 in Fig. |l|. The chiral phase transition is the second order for Nf = 2 in the chiral limit, while 
it is the first order for Nf = 3 and for m < nic = 0.0265. The mesonic masses can be defined as 
Q.{S) =Q.Q + ^M^o"^ + 2M2^„7r"2 _^ . . . ^jj^ S = ^+ X"{a" + in") /V2 parametrized with U(A^) 
generators (tr(A"A'') = 25"^). The temperature dependence of the masses are shown in Fig. |2[ For 
pseudo-scalar flavor-nonsinglet masses we find the Gell-Mann-Oakes-Renner relation 

{^+m)'^Mp^=mI?{^+m)^-m{qq), (4.2) 

if we identify + m as the pion decay constant Z^. On the other hand, the flavor singlet-masses 
have an additional contribution from the anomaly term as 

- AM^ , m2,o = + AM2 (4.3) 

with AMq = NfZ/{(j) + m)^. The would-be Nambu-Goldstone mode becomes massive due to the 
coupling to the Ua(1) interaction AMq, and the mass gap is related to the (replaced) quenched 
susceptibility f, similarly to the Witten-Veneziano formula. 

The topological susceptibihty is obtained as Xtop = d^^{S{d);d)/dd^\g=o with the saddle 
point solution S{d) = + /T]o(e)A*'/\/2 for small 6, and we find [0, §] 

1 1 1 

= - + — . (4.4) 

^top 

Here f is the modified susceptibility defined in Eq. ( 3.10| ). The fermion coupling screens Xtop via 
the contribution 

_ I?m{^+m) _ Ml^{^+mf 

Nf - Nf ■ ^^-'^ 

Noting that the Q dependence can be absorbed into the quark mass term as me^^/^f , one can show 
the axial Ward identity 

m ,_ . 

Ztop = -^ZpsO-^(w>- (4.6) 



6 



Chiral phase transition in a random matrix model with three flavors 



Takashi Sano 




Since the pseudo-scalar meson in the flavor singlet channel has nonzero mass ( |4.3[ ) because of the 
UA(l)-breaking term, the pseudo-scalar singlet susceptibility remains finite in the broken phase in 
the chiral limit. Thus for the small but nonzero quark mass m, the decrease of Xiop follows the 
chiral condensate {qq) ~ with increasing t, which is clearly observed in Fig. 3. 



5. Summary 

We have presented a chiral random matrix model where the determinant interaction is incor- 
porated by summing up the topological zero modes with the binomial distribution [^. The model 
has a stable ground state unlike the model of [Q], and describes the Nf dependence of the chiral 
phase transition. At the same time it reproduces the physical temperature dependence of the topo- 
logical susceptibility. This model can be applied to the 2-1-1 flavor case at finite temperature and 
density [g]. 

This work is supported in part by Grants-in-Aid of MEXT, Japan (No. 19540269 and 19540273). 
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